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ABSTRACT 


A  random  variable  X  with  distribution  func¬ 
tion  F(x)  can  be  written  as  x  =  R  (a)  ,  where  a  = 
F ( x)  and  R  =  F  .  The  function  R  (u)  is  the 
(first)  representing  function  of  X.  For  certain 
selected  distributions,  this  representing  function 
can  be  easily  expressed  (e.g.  logistic,  Cauchy), 
though  in  general,  approximation  or  tabulation  is 
required  (e.g.  Gaussian,  slash). 

A  situation  {X^  :i=i,  ...,  n}  is  a  collection 
of  independently  distributed  random  variables.  If 
the  are  identically  distributed,  the  situation 
is  termed  single,  otherwise  situation  is  termc 

compound .  For  example, 

X.  -  (l-«)F(x.)  +  <G(x.) 

is  a  simple  situation,  whereas  for  <=k/n, 

k  =  C , 1  , . .  .  ,n 

( 1  )  n  X '  s  ~  F  (  x) 


X 1  s  —  G(x) 


is  a  compound  situation. 

For  simple  situations,  the  low-order  moments 
of  the  order  statistics  can  be  conveniently  com¬ 
puted  in  terms  of  the  (first)  representing  func¬ 
tion  of  X.  For  compound  situations,  a  first 
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repr esen t i r.g  function  is  not  sufficient  for  com¬ 
puting  these  moments*  This  paper  provides  a  con¬ 
venient  method  of  computing  the  low-order  moments 
of  compound  situation  order  statistics  based  on 
hiqher  order  representing  function.  ihe  explicit 
derivation  of  the  second  representing  function  is 
given.  The  one-w i Id-Ga u s s i a n  situation  is  used  to 
illustrate  the  method.  Tables  of  one-wild- 
Gaussian  order-stati  st  ic  moments  are  displayed  for 
selected  sample  sizes. 
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1 .  In trod uc t i on . 

Order  statistics  play  an  important  role  in  statistics. 
Many  useful  estimators  are  based  on  linear  combinations  of 
order  statistics  (or  selected  subsets  thereof).  Informal 
inferential  procedures  (such  as  probability  plottinq)  are 
also  based  on  order  statistics.  Of  particular  importance 
are  the  low  order  moments  of  these  quantities,  specifically 
the  means,  variances,  and  covariances.  Tables  of  these 
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moments  exist  for  marv  of  the  comnor.lv  used  sampling  situa¬ 
tions.  In  almost  all  cases,  these  situations  are  "single", 
cor  re  3  pond  ir.q  to  a  sample  of  ir.de  pendent  and  identically 
distributed  random  variables.  (for  an  outstanding  except  ion 
see  Cavid,  Kennedy,  and  Knight,  197  7.) 

In  this  paper  v;e  provide  a  method  of  comp  it  ir.q  low- 
orcer  moments  of  order  statistics  from  “compound"  situations 
of  the  form 

n- 1  X  *  s  -  F  (  x) 

one  X  -  G  (x)  . 

The  Method  uses  what  we  call  the  second  represent  ir.q  func¬ 
tion  of  X,  namely 

:rVu)  U-o 

where  K  (u)  =  (u)  is  the  first  representing  function  of  X 

for  the  simple  situation 

X.  -  *W  =  (1  -*)E(x.)+«G (X.) ,  1  =  1,  n  . 

We  illustrate  the  method  using  the  one-wild  Gaussian  com¬ 
pound  situation 


n-1  X's  -  ’fr(x)  =  Gau(C,l) 

one  X  -  $(x/10)  =  Gau(0,100). 

Tli  i  s  compound  situation  has  been  used  extensively  in  studies 
of  robust/resistant  estimates  of  location.  The  case  where 
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vj  (x)  =  $(x/3)  uas  used  earlier,  and  is  one  of  the  cases 

tabulated  bv  Cavid,  Kennedy,  and  Knight  (1577).  tables  of 
the  low-order  moments  of  the  corresponding  order  statistics 
are  long  overdue. 

Section  2  cescribes  moment  calculations  for  simple- 
situation  order  statistics  in  terns  of  the  first  represent¬ 
ing  function.  Section  3  describes  moment  calculations  for 
compound  situation  order  statistics  in  terms  of  the  second 
representing  function.  The  one-w i Id-Ga us s ian  compound 
situation  is  used  to  illustrate  the  method  in  Section  4. 

2.  Simple  Situations. 

Consider  an  i id  sample  {x.:i=l,...,n}  of  random  vari¬ 
ables  with  distribution  function  t(x).  Let  v^  =  denote 

the  ith  order  statistic  with  v,  <  y~  <  ...  <  v  .  In  con- 
trast  to  the  x's,  the  v* s  are  neither  independent  nor  ident¬ 
ically  distributed.  Let  H(v.,Vj)  denote  the  joint  distribu¬ 
tion  function  of  v.  and  v ^ .  The  product  moment  of  v.  and  y^ 
>  y.,  is  given  by 

v . 

00  3 

■“ij  =  E<vivj)  =  I  J  vivj  dH(vi»vj) 

-  00  -  00 

where  dH(v.,v.)  is  proportional  to 

Fl“l(vl)£(vl)  [F(Vj)  -F  ( v.)  P”  1-1  f  ( v  j )  [1-F(v.)  ]  n"jdvldvj  .  (1) 

The  chanqe  of  variables  a  =  F ( v)  is  Monotone  so  that 
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\  -  ‘i-'ii  <  -  F'vj> 


me  the  above  expression  becomes 


u  . 
3 


m .  ■ 
ij 


=  f  T  R(u.)il(u.)dh(u.  ,1.) 

T  ’  i  J  i  ^ 

0  0 


where  dh(u.,Jj)  is  proportional  to 


ul*1  (Uj-J.)  i~  i_1  (1-Uj)  r,~^  du.du^ 


Thus,  given  the  representing  function  y  =  R(u),  the  low 
order  moments  can  be  obtained,  bv  numerically  integrating 

over  the  unit  triangle  0  <  u.  <  u.  <  1.  l\here  the 

—  I  —  j  — 

representing  function  cannot  be  qiven  explicitly,  a  numeri¬ 
cal  approximation  to  R(u)  is  required. 

*  a  special  form  * 

Quadrature  formulas  to  obtain  an  estimate  m.  ^  of  m.^ 
are  sometimes  more  convenient  if  the  region  of  integration 
is  the  unit  square  rather  than  the  unit  triangle,  and  if 
integration  involves  a  product  form  in  place  of  dli  (u.  ,’J ^  )  • 

This  is  easily  obtained  by  a  further  change  of  variables. 

Let 

u.  *  (l-z)w 

1  i 

1-Uj  *  (1-z)  (1-w) 

where  0  <  w  <  1,  c£z£1*  ?v^n  uj“Ui  *  2  ?rvl  since  th* 


diiSsuiaia 
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.  i  o n 


i-1  n  ,  r-i  j-i-1,,  .  n-  j+  i  .  ^ 

w  ( 1  — w )  JzJ  (1-2)  J  dwdz  . 

That  is,  w  and  z  are  independently  distributed  as  beta  ran¬ 
dom  variables: 


w  -  p(i,n-j+l)  and  z  -  p ( j- i , n- j- i+1 )  . 


The  alternate  expression  for 
is  therefore 


the  product  moment  of  v ^  and 


v . 
J 


1  1 

in..  *  f  J  a  (w(l-z)  )R  (v.(l-z) +z)  dpw(w)  d|E>2(  z)  . 

J  0  c 

If  desired,  one-d imen sional  quadrature  formulas  specialized 
for  integrating  a  function  of  a  beta- var  iable  could  now  be 
used,  iterating  the  integral.  The  accuracy  of  such  quadra¬ 
ture  rules  has  not  been  explored  in  detail. 


3.  Compound  Situations. 

Consider  a  realization  { x . : 1=1 , . . # ,n}  of  random  vari¬ 
ables  from 


n-k  X  *  s  -  F  (  x) 
k  X '  s  -  G  (x) 

for  kaG,...,n.  Let  v.  =  denote  the  ith  order  statistic 

with  v.  <  vn  <  ...  <  v  .  The  product  moment  of  v.  and  y.  is 
I—  2—  —  n  i  j 
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Blj  *  E‘vivj)  -  i  S  vivj 


-  00  oo 


u.iere  l;^(v.#Vj)  is  the  joint  distribution  of  v.  and  Vj  . 
This  joint  distribution  car.  be  derived  from  that  of 


VJ !  '  *  •  •  '  Vf.  ' 


!.(M(v.,v7,  2  k  J  n  3  (v.)  (r.-k)  1  n  F(v.) 

1  ;  n  UG  i^F 

k 


It  is  easv  to  see  that  the  resulting  formula  for  11^  (v^ 
is  appreciably  more  cumbersome  than  its  simple-situation 
counterpart*  Direct  integration  over  h^(v.,Vj)  is  not 
particularly  attractive,  especially  if  there  is  a  simple 
means  to  attain  the  sane  end* 


Consider  the  simple  mixture  situation 

Xi  ~  F<(xi)  =  (lH!t(xi)  +  «(xt)  1  =  l'  •••r  n  . 

The  joint  distribution  of  y,  and  y^  is  fc^(v. ,Vj)  and  can  be 
obtained  using  equation  (1).  This  leads  to  the  simple- 
mix  ture- si  tua  ti  on  order  statistic  moments: 

oo  vj 

=  /  J  vivj  dH<(vi'vj) 

-  oo  -  00 
1  uj 

*  J  J  R<('Ji)R<(ui)dK(ui,u.)  . 

0  C  J  J 

where  R^Cu)  =  F~M'J)  is  a  first  representing  function  for 
tbs  mixture. 
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V  . 

DO  3 


m  .  .  (<  ) 
13 


=  I  I  vivj  ^vV 


-  00  -  00 


n  oo  vj 

=  5  (^)  -<k { i  — € ) n_k  r  r  v . v .  dhk(v.,v.) 

k=n  ^  1  J  13 

K  u  -00-00 


-  5  (")  rij’  . 

k=0  K  ‘J 


This  fundamental  relationship  between  mixture  and  k-wild 
order  statistic  moments  allows  the  latter  to  be  calculated 
3  imply*  In  particular,  for  k  =  l,  equation  (2)  becomes 

n^H)  =  +  r.<  }  +  0  «2 )  . 

Differentiation  with  respect  to  <  and  evaluation  at  <  =  0 
lead  s  to 


^ij {<)  U-o  =  "n*ii) 


+  r.*n!P 


This  implies  that  the  one-wild  product  moment  can  be  written 
as  a  linear  combination  of  the  uncontam inated  product  moment 
and  a  term  due  to  the  contam  inat ion  viz 
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Algebraically,  the  correction  factor  is  obtained  bv  dif¬ 
ferent  iatinq  equation  (2): 


A*nij (<) 


1  V  > 


0  C 


<'  : 


arid  then  setting  <  =0 ,  to  give 


In  this  latter  equation  R,,  (u)  is  the  first  represent  ino 

u  ”  "  - 1 

function  of  X  at  <-G  contamination,  and  K  (u)  is  the  secor.c 

\  I 

representing  function  defined  bv  ^R^(u) 

Th  e  correspond  ing  conpo  und-  s  i  tuat  ion  o  rd  e  r-  st  a  t  i  $  t  ic 
moment  s  are  obtained. as 


( 1 J  , 
i  j 


1  a5 


-  ”S?1+k-,t  1  {Vw11  )+R  ,u  )R  (u j)'isMu1,-.1  ).  m 
0  0  L  J 


(_4ote  that  the  sample  size  enters  the  second  term  through 
both  n  and  dh(u^,Uj).)  This  expression  can  be  numerically 
evaluated  with  little  extra  effort  beyond  that  for  the 
simple-situation  moments  mj^.  Extensions  to  k-wild  com¬ 
pound  situations  are  easily  obtained  as  functions  of  the 
representing  functions  of  order  up  to  k+1,  where  in  general 
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4.  The  une-Mlc- jwVUisiar.  Situation. 

l\e  row  illustrate  tiie  preceding  discussion  using  trie 
compound  situation 

r- 1  X  *  s  -  $  (>0 

one  X  -  b  ( x/lC )  • 


To  do  so,  we  need  an  expression  for  the  first 
function  R^(u)  =  (u)  where 

a  =  ^(x)  =  (1— <)$(x)  +  <f>  ( x/1  0  )  . 


Now  sirice  R($(x))  =  x,  we  have 


K(^(x))  =  R(Mx))  +  (^(x)  )<sQ  +  0  (< 


=  X  +  -< *  r  (  f>  ( x )  )  *  f|>  ( x./l 0 ) - i>  ( x) j  +  0«2) 


where  r  ( u)  is  the  sparsity  function  —R  (u)  ;  see 
al  ( 1 S  4  7 )  for  the  original  definition.  For  our 
only  note  that  r  {  u)  is  easily  obtained  as 


_  Tdb  (xp^1 


r(u)  L  dx  Jx-R(u)  4(R  < u)  ) 


Hastings  et 
purposes  we 


2 

In  order  to  obtain  an  expression  for  R^(u)  as  x  +  0  (<*~)  , 
we  introduce  H[R(J>^(x))]  =  H(x) +(?(■<)  with 


Match  ^  , 


1CS1 


■  l(x)  =  -  r  (  Kx)  )  [  Kx./io  }  -  t>  (:<)  ]  . 

Tli  i  5  lead  3  to 

n(^(x)  )+<ll  [R  (^(x)  ))  =  x+OK2) 
or 

x  =  R  (u)  =  R(u)  +  <Ii  [R  (u)  ]  +  nj<2)  . 

The  first  and  second  re  pr  e  se  nt  i  nq  function^  of  X  are  nou 

% 

‘sas ilv  ootair.ed  as 
Rc(  u)  =  R  Co) 

H  ^  (  u )  ■Ii[il(u>]  =  -r(u)  U  (R  (u)  /lC)-u}  . 

j  (R  (u)  /10  )  -  i 
<5(R(u)) 

The  o r. e- wile  order  statistic  moments  can  now  be  numerically 
evaluated  bv  substituting  R ^ ( u )  and  (u)  into  equation  (3) 
Results  of  this  are  displayed  in  Table  1.  Xe  list  the  mean 
ana  covariances  of  the  one-wild  order  statistics  for  sample 
of  size  n  =  2  (1 ) 1 C •  For  comparison  purposes,  the  pure- 
Gaussian  order- sta ti st ic  moments  are  displayed  in  "able  2. 
As  expected,  the  effects  of  contain  ina  t  ion  are  most  strongly 
evidenced  in  the  extreme  (or  end)  order  statistics.  More 
detailed  tables  have  been  computed  by  A.  Bruce  (19CC). 
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